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ABSTRACT 

Wc present an analysis of the relationship between molecular gas and current star formation rate 
surface density at sub-kpc and kpc scales in a sample of 14 nearby star-forming galaxies. Measuring the 
relationship in the bright, high molecular gas surface density (Sh2 ^20 M© pc~^) regions of the disks 
to minimize the contribution from diffuse extended emission, we find an approximately linear relation 
between molecular gas and star formation rate surface density, Nmoi ~ 0.96 ± 0.16, with a molecular 
gas depletion time, ~ 2.30 ± 1.32 Gyr. We show that, in the molecular regions of our galaxies 

there are no clear correlations between r^"' and the free-fall and effective Jeans dynamical times 
throughout the sample. We do not find strong trends in the power-law index of the spatially resolved 

molecular gas star formation law or the molecular gas depletion time across the range of galactic stellar 
masses sampled (M, -lO'^ '^ - lO"-^ M©). There is a trend, however, in global measurements that is 
particularly marked for low mass galaxies. We suggest this trend is probably due to the low surface 
brightness CO J = 1 — 0, and it is likely associated with changes in CO-to-H2 conversion factor. 
Subject headings: galaxies: general — galaxies: spiral — galaxies: star formation — galaxies: ISM — 
ISM: molecules 



1. INTRODUCTION 

It is well known that the star formation activity in a 

galaxy correlates strongly with its stellar mass or stel- 
lar surface density (Dopita & Ryder 1994; Hunter ct al. 
1998; Kauffmann et al. 2003; Bhtz & Rosolowsky 2004, 
2006), and only weakly on galaxy morphology (Boselli 
et al. 2001). Characterizing the relations between stellar 
mass, star formation rate (SFR), and gas densities over 
cosmic time will provide important constraints on galaxy 
evolution by connecting the past history, present activ- 
ity, and future growth of a galaxy (Schiminovich et al. 
2010). 

A key factor in galaxy evolution is the rate at which 
molecular gas is converted to stars. Observational studies 
find that the relationship between SFR and gas content 
in galaxies can be written in the form, Ssfr — A Sg^^^, 
where Ssfr and Sgas are the SFR and gas surface densi- 
ties respectively, A is the normalization constant repre- 
senting the efBciency of the process, and N is the power- 
law index (Schmidt 1959; Kcnnicutt 1989, 1998). The 
gas can be atomic (HI) or molecular (H2) or a combi- 
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nation of both (HI-I-H2). This relationship is generally 
known as the Kennicutt-Schmidt law or the star forma- 
tion law, and variations on it are used as a empirical 
recipe in galaxy modeling (Schaye & Delia Vecchia 2008; 
Lagos et al. 2010). 

In recent years has become increasingly clear the im- 
portance of tmderstanding the link between star forma- 
tion activity and gravitationally bound molecular clouds 
(Wong & Blitz 2002; Tutukov 2006; Kennicutt et al. 
2007; Bigiel et al. 2008; Leroy et al. 2008, hereafter L08; 
Blanc et al. 2009; Verley et al. 2010; Onedara et al. 2010; 
Bigiel ct al. 2011; Schruba et al. 2011). Recent observa- 
tional studies indicate a strong relationship between the 
HI-to-H2 transition and the gravitational potential of the 
stellar disk and thus evince the connection between stel- 
lar pressure and the phase transition in the interstellar 
medium (Wong & Blitz 2002; Blitz & Rosolowsky 2004, 
2006). Numerical simulations that include cold gas and 
heating from young stellar population in evolution of in- 
terstellar medium (ISM) are in general agreement with 
the observations (Robertson & Kravtsov 2008). 

Most properties of galaxies, not least their star for- 
mation activity, are strongly correlated with mass (see 
Gavazzi 2009 for a review). Because new stars are 
born inside giant molecular clouds (GMCs), which them- 
selves evolve within the existing galactic stellar potential 
(Rafikov 2001; Li et al. 2005, 2006), it is important to 
characterize the interconnections between stellar mass, 
H2, and SFR (see Shi et al. 2011 for a recent study). To 
this end, we investigate the Eh2-5^sfr relation at sub-kpc 
and kpc scales in 14 nearby star-forming disk galaxies. 
While gas-SFR surface density relation provides an un- 
derstanding of on-going activities in the disk, the molec- 
ular gas depletion time {t^^) provides a measure of its 
future evolution. It is a quantitative measure of the effi- 
ciency of the star formation activity in molecular clouds, 
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defined as the time required for available molecular gas 
to be converted into stars while maintaining the existing 
rate of star formation (Roberts 1963; Larson et al. 1978; 
Kennicutt 1983; Kennicutt et al. 1994). 

In this study we measure in a sample of galaxies 
from the Survey Toward Infrared-bright Nearby Galax- 
ies (STING), observed in CO with the Combined Array 
for Research in Millimeter- wave Astronomy (CARMA), 
to investigate its relation with stellar mass and various 
dynamical timescales associated with gravitational insta- 
bility. The STING sample is composed of 23 northern 
{S > —20°), moderately inclined {i < 75°), high metallic- 
ity (12 + log(0/II) > 8.1) galaxies within 45 Mpc. These 
blue-sequence star-forming galaxies (Salim et al. 2007) 
have been selected to have uniform coverage in stellar 
mass (10*-i < M*/M0 < 10"-^), star formation activi- 
ties, and morphological types (A. D. Bolatto et al., 2011, 
in preparation). This is the second in a series of papers 
dedicated to exploiting the CO STING data set. In a pre- 
vious study we investigated the impact of methodology 
on the determination of the spatially resolved molecu- 
lar gas star formation law in NGC 4254, a member of 
the STING survey (Rahman et al. 2011; hereafter Pa- 
per I). In this study we explore relationships among cur- 
rent SFR, molecular gas, and stellar mass in normal star- 
forming galaxies using similar data analysis methodolo- 
gies as in Paper I. 

The organization of the paper is as follows. In |2] we 
briefly present multi- wavelength data set. In ^we pro- 
vide a description of the data products and data analysis 
methodology. Our main results and a general discussion 
are given i 4] and S|5] respectively. The conclusions are 
given in 36 

2. DATA 

The CO J = 1 — maps are obtained as part of the 
STING survey using the CARMA interferometer. While 
the full description of the CARMA observations will be 
presented in a forthcoming paper (A. D. Bolatto et al., 
2011, in preparation), we briefly mention here some as- 
pects of the interferometric data. The maps have a full 
sensitivity field-of-view of ~2' diameter, and are per- 
formed using a 19-pointing hexagonal mosaicing pattern 
with ~26" spacing. The angular resolution, measured 
as the full width at half maximum (FWHM) of the syn- 
thesized beam of the interferometer, varies from galaxy 
to galaxy in the range 6 ^3"-5". To construct molecular 
gas column density maps, we use the CO J ~ 1 — spec- 
tral cube to produce integrated CO J = 1 — intensity 
maps. To avoid integrating over noisy channels we use 
a velocity masking technique similar to that used by the 
BIMA Survey of Nearby Galaxies (SONG; Heifer et al. 
2003). 

To construct SFR tracer maps we use mid-infrared 
(MIR) 24 /im images from the Multi-band Imaging Pho- 
tometers (MIPS; Rieke et al. 2004) instrument on board 
the Spitzer Space Telescope. The calibrated images 
are obtained from the Spitzer Heritage Archive. We 
further process them by masking bright point sources 
such as Active Galactic Nuclei (AGN) which sometimes 
dominate the emission in the galaxy centers. The im- 
ages of NGC 1569, NGC 1637, NGC 3147, NGC 3949, 
NGC 5371, and NGC 6503 contain either foreground 
or background point sources, which we mask after vi- 



sual inspection. The central regions of NGC 1637 
and NGC 3198 contain bright sharply peaked point-like 
sources, despite the fact that galaxies are not known 
as AGN in the literature. The shapes of these cen- 
tral sources are consistent with the point spread func- 
tion of the MIPS instrument, which has a FWHM of 
^ 6". To remove the contributions of the point sources 
we mask the central regions of these two galaxies. Four 
STING galaxies (NGC 3147, NGC 3486, NGC 4151, and 
NGC 5371) harbor AGN at their centers, which we also 
masked. 

The maps of stellar mass are constructed from Two 
Micron AU Sky Survey (2MASS; Skrutskie et al. 2006) 
NIR 2.2 ^m(Ks-band) images. The image resolution is 2" 
for Kg-band images. At this wavelength masking of fore- 
ground and background objects was necessary for several 
STING sources. 

3. DATA ANALYSIS 

We select 14 galaxies from the 23 star-forming galax- 
ies in the STING sample. The selection is based on 
secured CO J = 1 — detection, availability of data at 
other wavelengths, and a H2 surface density threshold 
as we explain below. The subset spans stellar masses in 
the range 10^'^ < M^/M© < 10^^'^. Basic information is 
provided in Table [T] 

3.1. Molecular Gas and Star Formation Data 

We carry out our data analysis at two different reso- 
lutions: a fixed angular resolution (6"), determined by 
our SFR indicator, and a fixed spatial scale (1 kpc), 
to remove biases introduced by the different distances 
to our objects. The high resolution near-infrared and 
CO images were Gaussian-convolved to have the same 
angular as the 24 /j,m images. The 6" angular reso- 
lution covers a range of physical scales in our sample, 
TZ ^ 160 — 1250 pc, corresponding to physical distances 
D - 5.5 - 43.1 Mpc. The galaxies NGC 772, NGC 3147, 
NGC 4273, and NGC 5371, shown in bold face in Table 
[1] have distances such that 6" corresponds to 7?. > 1 kpc 
m spatial scale — those galaxies are left at their native 
resolution in the 1 kpc analysis, while the remaining 10 
galaxies are Gaussian convolved to the angular resolution 
corresponding to 1 kpc. All the analysis is carried out in 
images that have been regrid to the same pixel-scale and 
sampled at Nyquist-rate. We note here that, although 
we approximate the PSF of the MIR map as a Gaussian, 
the actual shape of the PSF is complex. It has prominent 
first and second Airy rings, with the second ring stretch- 
ing out to ~20". Nevertheless, approximately 85% of 
the total source flux is contained within the central peak 
with the FWHM of 6" (Engelbracht et al. 2007). 

We construct molecular gas surface density (Sh2) maps 
and SFR surface density (S24 nm) maps in a similar man- 
ner as described in Paper I. in the latter case following 
the prescription by Calzetti et al. (2007) who show that 
MIR 24 /im emission can be used as a SFR tracer for 
galaxies of normal metallicity and where the energy out- 
put is dominated by recent star formation. The MIR 
24 /im SFR tracer is given by. 
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Table 1 

Basic Information for tiic STING Sample 



Object 


Hubble 


D 


D25 


n 


5 


M, 




= mol 
~dcp 


log A 


N,„oi 


log A 


N„oi 


Corr. 


Corr. 




Typo 


Mpc 


kpe 


pe 


kpe 


Me 


Me 


Gyr 


6" 


6" 


1 kpe 


1 kpe 


P 


P 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 



NGC 628 


SA(s)e 


7.30t 


21.23 


211 


4.24 


9.95 


7.81 


2.91±0.75 


-0.23±0.05 


0.86±0.04 






+0.31 


-0.24 


NGC 772 


SA(s)b 


32.00 


51.77 


925 


18.62 


11.39 


9.71 


7.20±1.59° 


-0.59±0.19-'' 


0.88±0.10° 


-0.59±0.19 


0.88±0.10 


-0.40 


-0.67 


NGC 1637 


SAB(rs)e 


10.92 


11.38 


316 


6.35 


9.83 


8.25 


2.62±0.84 


0.07±0.1l'' 


0.69±0.08'' 


-0.20±0.22 


1.16±0.13'= 


-0.54 


-0.66 


NGC 3147 


SA(rs)be 


43.13 


43.39 


1247 


25.09 


11.36 


9.90 


6.16±1.35° 


-0.51±0.10-'' 


0.87±0.05° 


-0.51±0.10 


0.87±0.05 


-0.25 


-0.62 


NGC 3198 


SB(rs)c 


13.70+ 


21.11 


396 


7.97 


10.25 


8.14 


1.18±0.28 


-0.34±0.16 


1.16±0.09 


0.28±0.03 


0.78±0.02 


+0.25 


+0.09 


NGC 3593 


SAO/a?(s) 


5.50 


5.12 


159 


3.20 


9.81 


8.21 


1.98±1.00 


-0.60±0.09 


1.20±0.04 


-0.65±0.03 


1.23±0.02 


+0.37 


+0.14 


NGC 3949 


SA(s)bc? 


18.41 


11.74 


532 


10.71 


10.00 


8.22 


5.16±1.61 










+0.34 


-0.25 


NGC 4254 


SA(s)e 


16. 60+ 


24.33 


480 


9.66 


10.48 


9.44 


3.71±0.94° 


-0.24±0.06-'' 


0.86±0.03° 


-0.20±0.11 


0.82±0.05 


-0.50 


-0.61 


NGC 4273 


SB(s)e 


33.68 


18.50 


974 


19.59 


10.35 


9.34 


1.41±0.36 


-0.01±0.06 


0.93±0.03 


-0.01±0.06 


0.93±0.03 


-0.24 


-0.34 


NGC 4536 


SAB(rs)be 


14.40+ 


20.66 


416 


8.38 


10.24 


8.95 


1.13±0.37 


0.04±0.10 


0.96±0.05 


-0.03±0.11 


0.98±0.06 


-0.20 


-0.25 


NGC 4654 


SAB(rs)ed 


14.10* 


17.55 


472 


8.20 


10.21 


8.77 


2.11±0.54 


-0.19±0.05 


0.93±0.03 


-0.18±0.07 


0.84±0.04 


-0.16 


-0.37 


NGC 5371 


SAB(rs)be 


35.25 


39.89 


1019 


20.51 


11.30 


8.94 


4.34±1.16 


-0.25±0.09 


0.75±0.06 


-0.25±0.09 


0.75±0.06 


-0.61 


-0.65 


NGC 5713SAB(rs)be pee 


29.40+ 


22.63 


850 


17.10 


10.59 


9.59 


1.67±0.85 


-0.32±0.06 


1.09±0.03 


-0.29±0.11 


1.05±0.05 


+0.15 


0.00 


NGC 6951 


SAB(rs)be 


22.83 


23.53 


660 


13.28 


10.77 


9.37 


3.04±1.33 


-0.43±0.06 


1.00±0.03 


-0.50±0.07 


1.03±0.03 


+0.10 


-0.17 



Note. — Column (1): Galaxy name. Column (2): Galaxy Morphology from the RC3 catalog (de Vaucoulcurs et al. 1991). Column (3): 
Distance (D) in Mpc (^Prcscott et al. 2007; "^Tully et al. 2009). Column (4): Optical diameter (D25) in kpe. Column (5): Spatial resolution (7^) 
in parsec corresponding to 6 ' angular resolution at the adopted distance. Column (6): Physical extent (5) of the central arc-minute (in diameter) 
of the disk in kpe. Column (7): Logarithmic stellar mass (Mq). Column (8): Logarithmic molecular gas mass in Mq estimated within S including 
the contribution of helium. Column (9): Average molecular gas depletion time (t^^) and its Icr dispersion in Gyr at 6" resolution. Column 
(10) (11): Normalization constant and power-law index (Nmoi) of molecular gas SF law at 6" resolution. Column (12) & (13): Normalization 
constant and power-law index at 1 kpe resolution. Column (14): Spearman's rank correlation coefficient (p) from - tJ^°^ relation at 6" 

resolution. Column (15): Spearman's rank correlation coefficient {p) from relation at 6" resolution. 

The 6'' resolution of galaxies NGC 772, NGC 3147, NGC 4273, and NGC 5371, shown in bold face, corresponds to 7^ > 1 kpe in spatial scale. 
Stellar mass (M*) in Column (7) is derived from Two Micron All Sky Survey (2MASS) Kg-band total magnitude with a range of mass-to-light 
ratio, Tf^ ~ (0.3 - 0.8) M© /L 0,Ks, varying from galaxy to galaxy. It is this M* what we eall as the integrated stellar mass whieh incorporate 
the flux from the entire optical disk. The uncertainty is few percent in and ■ 

Parameters are derived for the surface density limit, Sh2 ^ 70 Me pc~^. 
^ Parameters are derived from Eh2 and S24/im maps with masked central region. The corresponding parameters for the unmasked central region 
are log ~(-1.20±0.09) and N„,ol ~(1.62±0.05). 

Central region was not masked at 1 kpe resolution. 



where L24 (in ergs^^Hz^^) is the luminosity spec- 
tral density at 24 /xm. We use a conversion factor 
Xco — 2.0 X 10^° cm~^/(Kknis~-^)to determine H2 from 
CO J = 1 — 0. The sensitivity (Icr) of the Sh2 maps 
varies among galaxies from '-^l.O to 5.1 Mq pc~^. The 
surface densities are multiplied by a factor of 1.36 to ac- 
count for the mass contribution of helium. 

3.2. Stellar Masses and Surface Densities 

Construction of stellar surface density maps (S*) in- 
volves two steps. First, we convert Kg-band luminosity 
(LkJ to stellar mass (M*). We use the B — V opti- 
cal colors for our galaxies and the relations by Bell & 
de Jong (2003) to compute the mass-to-light ratios in 
each galaxy. In our sample, we find a mass-to-light ratio 
Tf=~ 0.3 - 0.8 Mq/Lk,. Next, each pixel of the Lk, 
map of a given galaxy is multiplied by the appropriate 
T^' to derive the stellar mass map, i.e., M* = T^'Lk^. 
Each pixel value of the map is then divided by the 
area of the pixel to obtain the map. The galaxy inte- 
grated stellar masses are determined by integrating their 
H-hand 2-D surface brightness profile after masking of 
foreground or background objects. The NIR and MIR 
images were background subtracted prior to analysis. All 
surface densities have been inclination-corrected by ap- 
plying a cos i factor. 

There are uncertainties associated with the measure- 
ments of stellar masses. Although NIR emission is a 
good tracer of old (Gyr) stellar populations and it expe- 
riences little internal extinction, there is uncertainty as- 
sociated with the determination of mass-to-light ratios. 



For one, extinction affects the optical colors used in the 
calibration and consequently the employed mass-to-light 
ratio. More fundamentally, the mass-to-light ratio de- 
pends on the star formation history of the disk. Hence it 
varies considerably among galaxies of the same Hubble 
type and even within galaxies. For example, young M- 
supergiants near the plane of the disk (Aoki et al. 1991) 
as well as massive OB-associations (Regan & Vogel 1994) 
can contribute about ^ 50% of the total NIR disk emis- 
sion (Rix & Rieke 1993; Regan & Vogel 1994), potentially 
lowering the actual mass-to-light ratio and causing us to 
overestimate masses by factors of < 2. 

3.3. Analysis 

The sensitivity and flux recovery of the interferometric 
map limits the physical extent of the disk that can be 
studied. To minimize the effects of deconvolution and 
flux recovery problems we analyze the central 1' (in di- 
ameter) of the disks. Table [l] shows the physical extent 
of disk (S) corresponding to the inner arcminute for the 
STING galaxies. Following the methodology developed 
in Paper I we set a surface density threshold of SH2 ^ 20 
Mq pc~^. Studying the emission from the bright regions 
of the molecular gas and SFR tracer maps ensures that: 
1) the signal-to-noise is good, 2) interferometric deconvo- 
lution issues are minimized, 3) the potential contribution 
by the diffuse emission (DE) is less problematic, and 4) 
we focus on regions dominated by molecular gas. The 
DE is a component of the total disk emission that is 
unrelated to star formation activity and extended over 
the disk, in comparison to the localized emission associ- 
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ated with star formation. For example, a potential con- 
tributor to the DE at 24 /xm is infrared cirrus emission. 
The CO J = 1 — distribution of a galaxy can also con- 
tain DE not necessarily associated with the star-forming 
molecular clouds (Magnani et al. 1985; Blitz & Stark 
1986; Polk et al. 1988). Below our chosen surface den- 
sity limit the DE in both SFR and molecular gas tracers 
has the potential to affect the molecular gas-SFR surface 
density relation (see Paper I for a detailed discussion on 
this issue). 

Isolating the contribution of emission related to star 
formation activity from widespread DE is a complex is- 
sue. Recent work has shown that the DE has a poten- 
tially important impact in the determination of the star 
formation law, depending on its magnitude relative to 
emission coming from the star formation activity (Pa- 
per I; Liu et al. 2011). In Paper I we used an unsharp 
masking technique to remove a diffuse extended compo- 
nent, and showed that the most robust measurements of 
the star formation law are those performed on the bright 
regions of the studied galaxy (NGC 4254). In those re- 
gions the contribution from DE is least significant, and 
the recovered star formation law was approximately lin- 
ear. Here we adopt the methodology of Paper I, and 
minimize the impact of DE by focusing on the bright re- 
gions of our sample of disks. Finally, observational stud- 
ies suggest that Hl-to-H2 phase transition occurs around 
10-15 Mq pc-2 (Wong & Blitz 2002; Bigiel et al. 
2008; Leroy et al. 2008) and the nature of total gas-SFR 
surface density relation changes dramatically around this 
range (Kennicutt et al. 2007; Bigiel et al. 2008; LOS; 
Schruba et al. 2011 ). By focusing on high molecular gas 
surface density regions we avoid this issue. 

We carry out a pixel-by-pixel analysis (where the pix- 
els are sampled at Nyquist-rate) to probe the highest 
possible spatial resolution. We use the Ordinary Least 
Square (OLS) bisector method (see Isobe et al. 1990) to 
fit the molecular gas-SFR surface density relation and 
account for the measurement errors in each variable. As 
discussed in Paper I, it is important to consider the im- 
pact of the sampling methodologies, fitting procedures, 
and measurement errors when determining the I]h2-5^sfr 
relation. Some authors normalize Eh2 at a characteristic 
surface density prior to regression analysis to minimize 
the covariance between the exponent and the normaliza- 
tion constant of the power-law (see Bigiel et al. 2008 and 
Blanc et al. 2009 in this regard). We do not apply such 
normalization, and we have verified that our results are 
robust to this choice. 

In this study we use the 24 /im emission as the SFR 
tracer because it has several advantages over other single- 
or multi-wavelength tracers employed in the literature 
(Kennicutt et al. 2007; L08; Blanc et al. 2009; Verley 
et al. 2010). Most importantly, 24 fj,m images of uni- 
form quality are available for every galaxy in our sample, 
and no internal extinction correction is needed for this 
tracer. This is a major drawback for other SFR tracers 
at shorter wavelengths. Second, among various tracers 
studied in Paper I we find that the SFR obtained from 
24 /Ltm displays the tightest correlation with the molecu- 
lar gas. The scatter varies in the range ~ 0.1 — 0.3 dex 
depending on the subtraction of diffuse emission. Third, 
there is a striking spatial correspondence between the 
24 /zm and CO J = 1 — maps (see also Relaiio & Ken- 



nicutt 2009), suggesting that 24 fim is a faithful tracer 
of young (few million years), embedded star formation. 

In the following discussion, we use N,„oi to denote the 
power-law index of the molecular gas star formation law, 
where the H2 gas includes the contribution from helium. 

4. RESULTS 

At 6" resolution our analysis includes all 14 galax- 
ies. NGC 628 drops out of the sample at 1 kpc res- 
olution, however, as its peak surface density falls be- 
low the selected Eh2 threshold after smoothing. We 
have ^-^2000 ('-^lOOO) Nyquist-sampled pixel measure- 
ments arising from 14(13) galaxies at 6" (1 kpc) res- 
olution. The molecular gas and star surface densities 
span a wide range, Sh2 ~ 20 — 1000 M0 pc~^ and 
SsFR ~ 4 — 570 Mo Gyr~^ pc~^. Our results are pre- 
sented in Figures [T|4] The first three figures highlight 
spatially resolved cases whereas Fig. |4] shows global 
quantities such as average depletion time (r^p , mean 

value of Tj^^p' within the central arcminute, calculated in 
the logarithm), and integrated stellar mass as given in 
Table ttl 

4.1. Molecular Gas Star Formation Law 

Figure [1] shows the molecular gas-SFR surface density 
relation at 6" (top panels) and 1 kpc resolution (bot- 
tom panels). The diagonal dashed lines represent con- 
stant molecular gas depletion time defined as, r^l'^ = 

Sh2/Ssfr, assuming zero recycling of the materials by 
massive stars into the ISM. The left panels show the 
scatter plots where each Nyquist-sampled pixel has equal 
weight in the distribution. NGC 772 (shown in violet) 
and NGC 3147 (magenta) create a feature in the distri- 
bution with a slower rise in S24/im with I]h2 than other 



galaxies. We discuss this further in ^4.2 

A simple correlation test shows that the points in 
this diagram are strongly correlated: at either reso- 
lution we find Spearman's rank correlation coefficient 
p ~ 0.7. The OLS bisector method yields a power-law in- 
dex Nmoi ^l.liO.l at either resolution, where the error is 
derived from bootstrapping. Despite the fact these mea- 
surements are coming from wide variety of galactic envi- 
ronments and galaxy properties, the ensemble of points 
yields an approximately linear star formation law. 

To right side panels of Figure [T] provide a view of the 
sample that is not biased by galaxy size. In these panels 
any given measurement is weighted by the inverse of the 
total number of measurements of the galaxy to which 
the point belongs. Since larger galaxies contribute more 
points to the ensemble, the weighting scheme removes 
this bias by giving equal weight to each galaxy (Bigiel et 
al. 2011). The contours enclose 99%, 75%, 50% and 25% 
of the distribution. We note here that in extragalactic 
studies the inverse of rj^"' is sometimes known as the 
star formation efficiency (Young & Scoville 1991; McKee 
& Ostriker 2007). 



4.2. Molecular Gas Depletion Time 
We show Tj|"p' as a function of the molecular surface 
density I]h2 in the top panels of Fig. [2] The verti- 
cal hatch on the left of the panel demarcates the region 
where SH2 <20 Mq pc~^, the SH2 threshold discussed 
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Figure 1. Molecular gas-SFR surface density relation at 6" (top panels) and 1 kpc resolution (bottom panels) for the CARMA STING 
galaxies. The left panels show scatter plots where each pixel has equal weight in the distribution. The feature in the distribution below 
logjQ [S24^m/M0 Gyr~^ pc~^] < 1.4 showing a slow rise in S24^m with Sh2 is due to NGC 772 (violet open star) and NGC 3147 
(magenta open circle). The right panel shows a smoothed two-dimensional distribution where each point is weighted by the inverse of 
the total number of points of the contributing galaxy (all galaxies are equally important in the distribution, irrespective of the number of 
points they contribute). The contours of the smoothed distribution enclose 99%, 75%, 50% and 25% of the total. The diagonal dashed 
lines represent constant molecular gas depletion time t^°', in years. 



in ®] The individual measurements in STING galaxies 
are shown by green points, with the binned medians and 
(Icr) dispersions in black where the measurements from 
NGC 772 and NGC 3147 are excluded from the bins (see 
below). A horizontal dotted line represents the Hubble 
time, THubbie ^^13.7 Gyr (Spergel et al. 2007). 

It is clear from these diagrams that the H2 depletion 
time has at most very weak dependence on Eh2- A small 
correlation coefficient {p ~ 0.2) suggests that t^'^ is 
mostly uniform across the disk The measurements from 
NGC 772 and NGC 3147 are responsible for the plume 
of points showing a slow rise in the depletion time with 
I]h2 for Eh2 < 50 M0 pc~^. Indeed the I]h2-Ssfr re- 
lations in these two galaxies show flatter than average 



power-law indices (Nmoi ~ 0.8). The reason for this is 
unclear, but one possibility is that the contribution of DE 
to their 24 fim luminosity is worse than for the rest of the 
sample. Indeed Nmoi is increasingly closer to 1 if we in- 
crease the SH2 threshold over 20 Mq pc^^. These galax- 
ies are two of the most distant, thus most massive and 
intrinsically luminous, objects in our sample. Excluding 
NGC 772 and NGC 3147 we find r™°i ~ 2.30 ± 1.32 Gyr 

and r^°^ ~ 2.59 ± 1.39 Gyr at 6" and 1 kpc reso- 
lution, respectively (where the numbers correspond to 
the median and la dispersion in the logarithm of the 
measurements). The depletion time increases slightly to 
r™°i ~ 2.92 ± 1.87 Gyr and r'g"' - 4.17 ± 2.00 Gyr at 
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Figure 2. Molecular gas depletion time (t^°') versus molecular gas, Sh2i and stellar, S*, surface density at 6" and 1 kpc resolution. 
The horizontal dotted line represents the Hubble time, and the filled circles and associated error bars in black represent the median and 
la dispersion in Sh2 (S,) bins. Top panels: Results for Sh2- The gray hatch illustrates the region where Eh2 < 20 TsAq pc~^, which we 
remove from the analysis according to the discussion in jjs] NGC 772 (violet open star) and NGC 3147 (magenta open circle) are responsible 
for the plume of points showing a slow rise in the depletion time with Sh2 for Sh2 < 50 Mq pc~^. Bottom panels: Results for E«. The 
gray hatch corresponds to the region where E» < 80 Mq pc~'^, where t^°' increases for decreasing S, in the LOS study. NGC 772 and 
NGC 3147 are colored similarly in these panels as in Fig. ^ 



the respective resolutions when these two galaxies are 
included. In any case with the adopted definition t^'^ 
appears significantly shorter than the Hubble time. 

Our measurements agree very well with the results by 
Bigiel et al. (2011), who studied the molecular gas de- 
pletion time in 30 nearby (< 30 Mpc) spirals from the 
HERACLES Survey (Leroy et al. 2009). Given the min- 
imal overlap between the samples, the difference in the 
criteria used to select the galaxies in both surveys, and 
the fact that HERACLES is a single-dish CO J = 2 - 1 
survey, this agreement shows that r^"' determined in the 
molecular gas is very similar in disks across a wide range 



of galaxy properties. 

We show Tj^p' versus the stellar surface density in 
the bottom panels of Fig. [2j The hatched section on the 
left shows the region where t^°' decreases monotonically 
with increasing E*, and HI dominates the gas surface 
density (LOS). We do not find any correlation between 
'''dcp' ^^"^ E* showing that the molecular gas consumption 
time is independent of stellar surface density in the re- 
gion of the disk where H2 is the dominant component of 
the ISM (see also LOS for a similar result). Since stellar 
mass surface density dominates the underlying gravita- 
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Figure 3. Molecular gas depletion time (1"^°') versus local dynamical 



timescales. Left panel: t^°' versus Jeans time (t™"'). Right panel: 



r',"" versus free-fall time {t^°^). Results are computed at 6" resolution. The measurements of NGC 772 and NGC 3147 have similar color 



codes and legends as in Fig. jlj 
[ati< 

entire sample (see discussion in i]4.3| 



There is no correlation between t^°' and 



r™° , either on average or galaxy by galaxy. By contrast, there 
is a reasonably strong correlation between t^^°' and t™°' on a galaxy-by-galaxy basis that is somewhat washed out when considering the 



tional potential in these objects, this result suggests that 
locally the molecular gas-SFR surface density relation is 
independent of the large scale galactic potential. In other 
words, once diffuse molecular clouds turns into isolated, 
self-gravitating objects, the conversion of II2 into stars 
in GMCs is not sensitive to the overall gravitational po- 
tential GMCs reside. 

Note that while both star formation rate and molec- 
ular gas mass surface densities correlate with S*, their 
ratio T^l'^ does not correlate with either quantity. Both 

panels of Fig. [2I show r^""' is independent of I]h2 and 
E*. For the STING data set we find that the correlation 
coefficients for S24^m — E* and I]h2 — E* relations are 
p ~ 0.48 and p ^ 0.63, respectively, when all galaxies 
at 6" resolution are included. The correlation strength- 
ens, p ^ 0.65 and p ^ 0.70 for the respective relations, 
if the contributions from NGC 772 and NGC 3147 are 
removed from the distribution of points. The correlation 
coefficients are similar at 1 kpc resolution. 



4.3. Dynamical vs. Star Formation Timescale in the 
Molecular Disk 

Molecular gas and stellar mass surface densities as 
shown in Fig. [2j along with the gas and stellar velocity 
dispersions, can be used to derive dynamical timescales 
associated with the growth of GMCs such as the Jeans 
time (t™°0 and the free-fall {Tg°^) time. We would ex- 
pect SFR oc Sh2 /rf°' or I]H2/Tf °' if either of these dy- 
namical timescale is relevant to star formation, which 

.,_mol 

'dcp 



translates into a proportionality between t™°' and either 
dynamical time (LOS; Wong 2009). 

For a plane parallel, axisymmetric, and isothermal two- 
component (gas and star) disk under hydrostatic equilib- 



rium, the Jeans time can be written as. 



,_mol 



1 



Eh2 



1 



^H2 



(2) 



under the assumption that the gas surface density is 
dominated by II2. The quantities Cg and c* are the ve- 
locity dispersion of gas and stars respectively along the 
2-direction, and h^ g is the vertical scale height of the 
gaseous disk. The Jeans time defined this way should 
be regarded as the "effective" Jeans time since originally 
it was defined in terms of thermal motions, with Cg cor- 
responding to the sound speed. Our t™°' is defined in 
terms of the gas velocity dispersion, which is dominated 
by turbulent motions. The free-fall time of gas in the 
same disk can be written as, 



37r 1 
32 G^ 

73 

4G ^Eh2 



-'H2 



-1/2 



(3) 



where pH2 is the mid-plane H2 volume density, related 
to the surface density by pH2 = EH2/2hz^g. These rela- 
tions stem from the fact that in the two-component disk 
the scale height and the velocity dispersion of each com- 
ponent are interrelated (Kellman 1972; Talbot & Arnett 
1975; van der Kruit 1983, 1988). For example, the rela- 
tion between stellar velocity dispersion and scale height 
can be written as (see equation 16 in Talbot & Arnett 
1975), 



c* = ttG hz^* (E*/c* -I- I]H2/Cg). 



(4) 



In order to estimate dynamical timescales we require 
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Figure 4. The molecular star formation law as a function of galactic stellar mass. Top panel: Power-law index (N„joi) versus stellar mass 
at 6" resolution. The black points and associated error bars correspond to the determinations in our sample of galaxies. The horizontal 
dotted lines represent indices for linear (Nn^oi = 1) and non-linear (Nnioi = 1-5) molecular gas star formation laws. The panel also shows 
using a gray square ^rnol for NGC 1637 if its nuclear region i s no t masked. The measurements of NGC 772, NGC 3147, and NGC 4254 
shown in gray and black employ different Sh2 thresholds (see i |4.4| for details). Bottom panel: Molecular gas depletion time versus galactic 
(total) stellar mass. The average depletion time, "r^""', is calculated as the average of the logarithm of Sh2/Sspr in regions where Sh2 
is over the adopted threshold. The horizontal dotted line represents the Hubble time. The gray stars and circles show the global t^°' 
measurements within the optical disk of star forming galaxies taken from LOS and Kennicutt (1998), respectively. NGC 3949 has been 
shifted by -1-0.02 dex horizontally for clarity of presentation. 



both Cg and c* , but direct measurement of either quan- 
tity is challenging. We use the fact that c* and h^^ are 
interrelated and so to obtain c* we first estimate stellar 
vertical scale height using the empirical relation between 
radial scale-length (hR^*) and vertical scale-height (h^^*) 
fromKregelet al. (2002), < hR,*/h^,* >= 7.3±2.2, where 
hR„* for STING galaxies is olDtained from the 2MASS 
data. For Cg we adopt a constant value of 10 km sec~^, 
very similar to that used by LOB and Wong (2009). 
Figure Is] shows the correlations between r^""', t™"', 



and r^°' at 6" resolution. When the entire STING sam- 
ple is considered we find no correlation between t^'^ and 
either of the dynamical timescales (the correlation coef- 
ficients are p < 0.1 for either timescale). Wong (2009) 
analyzed several star-forming galaxies employing radial 
profiles, and finding weak or no correlation between these 
quantities. Our results, with slightly different normaliza- 
tion for T™°' and t^°' and with a general solution of c* , 
corroborate these results in the molecular regions of disks 
for the ensemble of STING galaxies. 
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Figure 5. Molecular gas depletion time versus specific star formation rate sSFR (in yr~^) showing both resolved (6") and global 
measurements. The resolved measurements are shown in faint gray in the background. The averages within the central arcmin and the 
associated Icr errors are shown with filled circles and thick black lines respectively. The legends are similar to Fig. |4] The horizontal dotted 
line represents the Hubble time and the diagonal dashed line is the fit log[r^°'] = —0.724 log[sSFR] + 1.54, from Saintonge et al. (2011). 
There are 7 outliers at the top- left corner associated with the spatially resolved measurements of NGC 5371. 



When considering individual galaxies, however, a clear 
negative correlation emerges between and r^°' for 
most galaxies in the STING sample. Table 1 shows 
that this result is significant: the correlation coefficients 
for r^p^ vs. r™°' are small and scattered around 
pretty much symmetrically, suggesting that the growth 
of Jeans-type dynamical instabilities is not responsible 
for the regulation of star formation. By contrast, for 
^dcp vs.T^°' half the sample has -0.7 < p < -0.5 and 
only one galaxy has p }^ 0- Although the sign of the 
correlation is puzzling, taken at face value this would 
suggests that some type of gravitationally driven accre- 
tion of molecular gas onto GMCs may be related to the 



regulation of star formation in the inner disks, and that 
the degree to which this plays a role may be different 
from galaxy to galaxy. But there are key caveats that 
make us skeptical of such conclusion. 

We think that the observed correlation is likely not 
physical, but introduced by how we calculate t™°' and c.^ , 
and the fact that moderate (p ^ 0.5) to strong (p > 0.7) 
correlations exist between Eh2, S*, and Ssfr- Thus 
there is a coupling between the X and Y axes computed 
in this manner that is difficult to avoid with the current 
state of the observations. To explore further we conduct 
a Monte Carlo experiment and use the fact that S24/xmj 
Eh2, and E* are interrelated (see the appendix for a 
detailed account of the numerical experiment). We ap- 
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proximately reproduce the observed correlations, vs. 
5^24 /xm and vs. Eh2, assuming E* as the independent 
variable. Our numerical experiment, although crude, 
strongly suggest that the observed galaxy by galaxy cor- 
relations may be almost entirely attributed to this ef- 
fect. Additionally, there is considerable uncertainty in 
the degree of coupling between scale heights and velocity 
dispersions that complicates the interpretation. 

4.4. Global Correlations 

The top panel of Figure |4] shows the distribution of the 
molecular power-law index Nj^oi as a function of galaxy 
(total) stellar mass M* . We only display the results at 6" 
resolution, as the results for a fixed 1 kpc spatial scale are 
effectively identical (see Table 1). Most of the galaxies 
show an approximately linear relation between molecular 
gas and SFR surface densities, with little indication of 
a trend with galaxy stellar mass. The STING sample 
average is Nmoi ~ 0.96 ± 0.16. Note that the panel does 
not show the measurement of Nmoi for NGC 3949, which 
has a very small dynamic range (~0.3 dex) in both Eh2 
and EgpR yielding a correlation coefficient close to zero. 

The largest value of the power-law index is found for 
NGC 1637, Nnioi ^ 1-6, if we do not mask its nuclear 
region (gray square in Fig. |4). As mentioned in section 
fp] despite not being a bonande AGN the 24 fim image 
of NGC 1637 contains a point source at the center. This 
source makes the Eh2 — SgpR, relation very steep when 
included. On the other hand, when the central part of 
the is masked in order to remove the contribution of the 
point source the power-law index becomes significantly 
flatter (Nmoi ~ 0.7). This illustrates the potential effect 
of AGN or nuclear clusters on the star formation law 
determination, and some of the care that should be exer- 
cised to obtain unbiased results (see also Momose et al. 
2010). 

We find significantly flatter than unity power-law in- 
dices (Nnioi ^ 0.7-0.8) in three other galaxies: NGC 772, 
NGC 3147, and NGC 4254. The indices for these galax- 
ies are shown in gray in the top panel of Fig. |4j The 
reason for the low index in the power-law fits is that 
in these sources E24^ni is almost constant in the range 
Eh2 ^ 20 — 70 M0 pc^^. Consequently the distribu- 
tion of points in the Eh2 ~ Ssppj plane possesses fiat 
tails at the low surface density end. Although the ex- 
tended tail contains a small number of points (< 20 per 
galaxy), these points significantly impact the regression 
analysis, yielding flatter indices. By contrast, we ob- 
tain Nmoi ~ 1 in each case by raising the surface density 
threshold to Eh2 ^ 70 Mq pc~^, showing that the high 
surface brightness regions display a linear trend. 

The average molecular depletion time, T^p, in the 
STING galaxies is shown in the lower panel of Fig. |4] 
as a function of galaxy stellar mass. The average is 
computed in the logarithm over the pixels where Eh2 
is larger than the threshold. For the STING sample we 
find f^°i - 3.19± 1.91 Gyr with a range of 1.1 to 7.2 Gyr 
for individual galaxies (see Table 1). This sample aver- 
age is slightly longer than the previous mentioned 
results for the spatially resolved case, because of the "per 
galaxy" weighting. NGC 772 and NGC 3147 are again 
the contributing galaxies that to push f^^p' to the larger 
value. 



We observe a weak correlation between r^^"' and galaxy 
mass: a linear regression analysis yields logf^^^p — 

(0.50 ± 0.13) logM, -I- (4.27 ± 1.38). We should cau- 
tion, however, that this correlation is entirely driven 
by the high-mass galaxies — particularly NGC 772 and 
NGC 3147. The rank correlation coefficient between f^l'^ 
and total M^, is p ^0.50 and ^0.15 with and without 
these two galaxies. As discussed in the previous para- 
graphs, these are the galaxies where we also see a small 
Nnioi- Note that although the derivation of stel lar m ass 
from NIR light has uncertainties (discussed in ^3.2) the 
large dynamic range in M.^ lessens their impact on this 
study. 

In a recent study, Saintonge et al. (2011) report an 
interesting correlation between H2 depletion time and 
stellar mass in star-forming galaxies, in the stellar mass 
range 10^° < M^/M© < 10^^-^. Using unresolved 
molecular gas measurements and estimating SFRs using 
FUV-l-optical spectral energy distribution fitting, Sain- 
tonge et al. (2011) find a global disk-average r^"' ~1 Gyr 
in star-forming galaxies which increases by ~ 0.6 dex 
for galaxies with stellar masses from M* ~ 10^° M© to 
M„ ~ 3 X 10^^ Mq. By comparison the mean resolved 
depletion time of the STING sample is f^°i ~ 3.0 ±1.6 
Gyr, consistent within the uncertainties with the deple- 
tion time found by both L08 and Bigiel et al. (2011) who 
employ the same methodology. 

The difference between these two sets of results is due 
to the fact that the studies measure fundamentally dif- 
ferent quantities. The resolved r^"' measurement is car- 
ried out in regions where CO is detected, and it is thus 
a statement about the depletion time in GMCs (partic- 
ularly inner disk GMCs). The unresolved t^"' measure- 
ment includes emission from all regions inside the beam, 
particularly the outer disk, and it includes SFR tracer 
emission from regions that have little or no molecular 
gas. It is thus a statement about the time it would take 
for the galaxy to run out of molecular gas, modified by 
variations in the Xco factor (which are likely to be sig- 
nificant for outer disks) and contribution from diffuse 
emission and regions of star formation that have no de- 
tectable CO emission to the global SFR (which can be 
very significant in low surface brightness regions). 

To explore further whether a f^p — M* correlation is 
a general property of star-forming disk galaxies we com- 
pile the global measurements from L08 and Kennicutt 
(1998). This sample extends the dynamic range in stel- 
lar mass by more than an order of magnitude compared 
to Saintonge et al. (2011), incorporating many more low 
mass systems 10^-^ < M^/M© < 10"-^ The measure- 
ments of L08 and Kennicutt (1998) are shown by open 
stars and open circles, respectively, in Fig. |4 L08 use 
Tf = = 0.5 M©/Lk, to derive the stellar mass tor a sam- 
ple of 23 galaxies. This value is roughly in the middle 
of the range of T^" used in this study. The SFR es- 
timation and Xco used by L08 are similar to ours, so 
no adjustment of the measurements are necessary. The 
panel shows the measurements of 19 galaxies in L08 that 
have secure Mhj estimates. Kennicutt (1998) provides 
the global Ssfr and Eh2 of a sample of 61 star-forming 
galaxies. We adjust the original SgpR and Eh2 measure- 
ments to be consistent those of L08 and this study. To 
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derive for these galaxies we use Tf^= — 0.5 Mq/Lk, ■ 
There is a minor overlap between the STING, and the 
LOS or Kennicutt samples. On the other hand, there is a 
substantial overlap between the samples of LOS and Ken- 
nicutt. Since we do not intend to make any quantitative 
relation between f^^^p and M*, we treat these samples 
as independent because of their significant methodolog- 
ical differences. It is apparent from Fig. |4] that there 
is a trend in the global measurements, once the small 
mass galaxies are included. The global t^°^ varies al- 
most three orders of magnitude in the stellar mass range 
10^-^ < M» < 10^^-^, where small mass galaxies have 
systematically lower global t^°^. The trend for the small 
mass galaxies is simple to understand in terms of a higher 
value for Xco in galaxies that are of lower mass and lower 
metallicity (Leroy et al. 2011; Krumholz et al. 2011). A 
recent study of the Small Magellanic Cloud that avoids 
using CO to trace II2 at low metallicity finds a strong 
metallicity effect in the CO emission, but not a measur- 
able one in r^"', suggesting that the trend in Fig. jZjis 
mostly an Xco effect (Bolatto et al. 2011). 

Are systematic changes in Xco associated with metal- 
licity strong enough to explain the magnitude of the ob- 
served trend? Reproducing the mean behavior of the 
unresolved data requires typically increasing Xco by a 
factor of 1.3 dex (^^20) going from galaxies with stellar 
mass log M* ~ 10.5 to 9.0. The corresponding metallicity 
change expected, according to the mass-metallicity rela- 
tion (Tremonti et al. 2004; Mannucci et al. 2010) is ap- 
proximately 0.6 dex. That is approximately the change 
in metallicity between the Milky Way (Baumgartner & 
Mushotzky 2006) and the Small Magellanic Cloud (Pagel 
2003). Our best estimates of Xco on the large scales in 
the Small Magellanic Cloud show that it is 30 to SO times 
larger than Galactic (Leroy et al. 2011). If this behavior 
is typical, we conclude that Xco changes could indeed be 
the main driver behind the trend for global t^""' apparent 
in the lower panel of Fig. |4j 

The past evolution of a galaxy is usually quantified by 
the the stellar mass assembly time (Kennicutt et al. 1994; 
Salim et al. 2007). It is the time required for a galaxy 
to assemble the current stellar mass at its present SFR. 
Assuming zero recycling of the materials by massive stars 
into the ISM, this timescale is defined as, — E^/Ssfr- 
The inverse of this timescale is commonly known as the 
specific star formation rate (sSFR). 

Figure [51 shows t^°' versus sSFR for the STING sam- 
ple as well as the literature data discussed earlier. It 
also shows the relation between and sSFR for 

star-forming galaxies in the local universe reported by 
Saintonge et al. (2011). Note that there are 7 out- 
liers at the top-left corner associated with the spatially 
resolved measurements of NGC 5371 which are com- 
ing from its central regions where each of pixel has 
>2000 M0 pc^^. Local galaxies form a bi-modal dis- 
tribution in the M* - sSFR plane (Salim et al. 2007) 
where star-forming galaxies form an horizontal branch. 
In terms of Fig. [5j this observation, however, would im- 
ply that the large, massive galaxies would form a locus 
at the upper left whereas the low mass system would 
fall to the lower right, likely due the Xco effects dis- 
cussed above. The global measurements from the liter- 



ature are broadly consistent with the fit by Saintonge 
et al. (2011), with significant deviations at low masses. 
The resolved 6" (gray dots), 1 kpc (not plotted), and 
average (black symbols and error bars) measurements of 
STING galaxies show a similar correlation between t^°^ 
and sSFR (p > —0.6). The differences in the normal- 
ization between these measurements and the empirical 
relation reflects the difference in the aperture selection 
e.g., central arc-min vs. extended disk. As it can be sur- 
mised from the fact that r^p is approximately constant 
with S* (Fig. [2]), the negative correlation observed in 
the STING data set is likely attributable to the fact that 
the horizontal axis is proportional to the SFR whereas 
the vertical axis is proportional to the inverse of the SFR. 
We have further explored this correlations using a Monte 
Carlo experiment which suggests that the negative cor- 
relation between t'^°^ and sSFR is governed at the fun- 
damental level by the inter-connections among E.^, Eh2, 
and E24/im (see appendix for more on this experiment). 

5. DISCUSSION 

In this section we interpret the results of this study 
and relate them to the existing scenarios of star forma- 
tion and the molecular ISM. In the spatially resolved 
case our results can be summarized as follows: within 
the range of Sh27 S*, and M* explored in this study, 1) 
the resolved molecular gas depletion time is independent 
of cloud properties in the disk such as Eh2- 2) Dynami- 
cal timescales, such as the (effective) Jeans time and the 
free-fall time in the molecular disk, do not correlate with 
the molecular gas depletion time over the entire sam- 
ple. And, 3) the resolved molecular gas depletion time 
is approximately independent of the disk environment 
represented, for example, by S* or M*. 

The uniformity of t^°^ over a wide range of Sh2 is 
most naturally explained as the consequence of the ap- 
proximate constancy of the depletion time for molecular 
gas in GMCs. Indeed, observations of galaxies in the 
Local Group and beyond suggest that the properties of 
GMCs are fairly uniform (Blitz et al. 2007; Bolatto et al. 
2008; Fukui & Kawamura 2010; Bigiel et al. 2010). In 
this scenario a linear star formation law follows naturally, 
where the Sh2 — Ssfr relation arises from the number 
of GMCs fiUing the beam (Komugi et al. 2005; Bigiel et 
al. 200s). A linear molecular gas star formation law is 
consistent with the scenario in which GMCs turn their 
masses into stars at an approximately constant rate, ir- 
respective of their environmental parameters (Krumholz 
& McKee 2005). 

For our second result we determined two of the dynam- 
ical timescales associated with the gravitational growth 
of GMCs in the molecular part of the disk; the effective 
Jeans time and the free-fall time. This is a challeng- 
ing determination, because it relies on a number of as- 
sumptions to estimate the relevant velocity dispersions 
and scale heights. The lack of a clear correlation be- 
tween these dynamical timescales and t™°^ suggest that 
the star formation regulation is not necessarily associ- 
ated with the growth time of large scale gravitational 
instabilities that may create and grow GMCs. We see a 
correlation between r^""' and t^°^ in individual galaxies. 
Indeed, approximately half the sample shows a correla- 
tion, although it is a physically puzzling inverse corre- 
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lation where longer free-fall times correspond to shorter 
molecular depletion times, thus more star formation ac- 
tivity per unit molecular mass. B\it our numerical exper- 
iments suggest that where we see it, it is explained by the 
observed correlations between Sh2) ^sfr, and E*. 

The third result connects t™"^ with the local gravita- 
tional potential traced by S*. This result implies that 
T^p is mostly independent of the local potential in the 
molecule rich regions of the disk. LOS reaches a simi- 
lar conclusion, and our result corroborates their findings. 
Following theoretical studies of Elmegreen (1989. 1993) 
and Elmegreen & Padoan (1994) that rely on the; equi- 
librium balance between H2 formation and radiative dis- 
sociation, Wong & Blitz (2002) and Blitz & Rosolowsky 
(2004, 2006) suggested that the mid-plane hydrostatic 
pressure plays a critical role in governing the equilibrium 
fraction of molecular gas phase in the disk. A recent the- 
oretical study by Ostriker, McKee, & Leroy (2010) ex- 
plains the same observations as due to the equilibrium in 
a multiphase ISM, where the stellar potential plays also 
an important role. If the disk gravity is dominated by 
the stellar potential and the gas scale height is smaller 
than the stellar scale height, these studies show that the 
molecular ratio, defined as i?nioi = 5]h2/5^hIi increases 
approximately linearly with the ambient pressure, 

R^oi a P« a (Sgas ^*°-'r, (5) 

where a ^ 0.8 — 0.9 and Sgas is the mid-plane total gas 
(HI-I-H2) surface density (LOS). In the regions that are 
dominantly molecular, however, i?moi ^ 1 and Sh2 ~ 
Sgas, thus Tj^°' is independent of E*. In other words, the 
stellar potential plays a role at determining the fraction 
of gas that is molecular, but it does not affect the rate 
at which GMCs collapse and their gas is converted into 
stars. Blitz & Rosolowsky (2004, 2006) show that the 
ISM becomes molecular around w 120 ±50 Mq pc~^. 
The mean (median) of the STING galaxies is 580 
(395) M0 pc~^. It is, therefore, very likely that we probe 
the molecular ISM where the assumption i?nioi 1 is 
well justified (Xue et al. 2011, in preparation). 

We should, however, stress here that whether the for- 
mation of molecular clouds is regulated by the mid-plane 
pressure (Blitz & Rosolowsky 2004, 2006), gravitational 
instability (e.g., Mac Low & Glover 2010), or photo- 
dissociation (Krumholz et al. 2009) is still not settled. 
A combination of dynamical and thermodynamic factors 
may be required to regulate the formation GMCs and 
the subsequent star formation inside the clouds (e.g., Os- 
triker et al. 2010). 

We find that the resolved molecular gas depletion time, 
averaged over the central regions of our galaxies, shows a 
positive but weak correlation with the integrated stellar 
mass in the stellar mass range lO'''^ < M*/M0 < 10^^ '''. 
Given that the correlation is dominated by a few galaxies 
in this sample, it is difficult to assert it with any degree of 
confidence. If it were real, it would suggest that there are 
weak environmental effects on the in the GMCs, and 
analytical models must be able to reproduce this large 
scale behavior. A systematic resolved study with a large, 
well-defined sample having homogeneous measurements 
over the extended disk is necessary to shed more light on 
this issue. 



Finally, using literature data wc show that star forming 
galaxies, spanning a large dynamic range in stellar mass 
(10^-^ < M*/M0 < 10^^-^), show a clear correlation be- 
tween molecular depletion time and galaxy mass. This 
correlation is, at least in the low-mass galaxies, most 
likely explained as arising from systematic Xco varia- 
tions, although systematic trends in the SFR calibra- 
tion may also play a role. Fundamentally, however, it 
highlights different physical processes than the resolved 
molecular gas measurements. The approximate con- 
stancy of r^p in the resolved molecular measurements is 
a statement about the depletion time in GMCs, which are 
regions of the disk dominated by molecular gas. The un- 
resolved measurement provides information about 
the time it would take for the entire galaxy to run out of 
molecular gas at its current SFR, and it folds in the ef- 
fects of the Xco factor and the contribution to the global 
SFR from emission arising in regions with no CO. 

6. SUMMARY 

We present a comprehensive analysis of the relation 
between Eh2 and Esfr as a function of stellar mass at 
sub-kpc and kpc scales using a sample of 14 nearby star- 
forming galaxies observed by the CARMA interferome- 
ter spanning the range of galaxy stellar masses 10^ '' < 
M»/Mq < 10^^'^. We measure the relation in the bright, 
high molecular gas surface density (>20 Mq pc^^) re- 
gions of star-forming disks. Sampling these CO-bright 
regions has the advantage of minimizing the contribu- 
tion from diffuse extended emission present in the SFR 
tracer and molecular gas disk. Our main results are: 

1. The pcr-galaxy average star formation law for the 
sample, determined using E24^m as the SFR in- 
dicator, is f^°i ~ 3.19 ± 1.91 Gyr and N^oi ~ 
0.96 ±0.16. 

2. The resolved molecular depletion time t^°' is inde- 
pendent of both molecular and stellar surface den- 
sities, Eh2 and E*, respectively. We find r^"°' ~ 

2.30 ± 1.32 Gyr and r^°^ ~ 2.59 ± 1.39 Gyr at 6" 
and 1 kpc resolution, respectively. 

3. There is no clear correlation between t^'^ and the 
effective Jeans time, rf°i, or the free-fall time, 
Tg°^, in the molecular regions of our galaxies. 
These dynamical timescales, which may be impor- 
tant for GMC growth, do not appear to regulate 
the star formation once the gas is molecular. 

4. There are no strong trends across our range of stel- 
lar masses for either the power-law index or the 
normalization of the resolved molecular star for- 
mation law. 
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APPENDIX 

A. INVESTIGATING CORRELATIONS BETWEEN MOLECULAR GAS DEPLETION TIME VS. DYNAMICAL TIMESCALES 

AND SPECIFIC STAR FORMATION RATE 

In this study we show that the observed trends in molecular gas depletion time vs. dynamical timescales for individual 
STING galaxies vary significantly, showing little or no correlation to strong negative correlation (see Table IT]) . The 
correlations for either dynamical timescales, however, vanishes (p ~ —0.01 in t^'^ - t™°' and p ^ —0.16 in t^'^- t^°') 
when the entire ensemble is considered (see Fig. Isl). The molecular gas depletion time, on the other hand, shows strong 
negative correlation (p ~ —0.64) with the specific star formation rate (see Fig. [5]). In this appendix we demonstrate 
by a simple Monte Carlo experiment that these trends, particularly t^'^ - sSFR relation, are mostly driven by the 
existing correlations among Eh2, ^2ifim, and and their connections to the stellar velocity dispersion. We carry 
out this experiment using the 6" resolution STING data set. 

We begin with the assumption that the stellar mass surface density is a fundamental parameter that influence 
both local and global evolution and organization of molecular gas in the disk and star formation. Indeed observational 
evidence shows that almost all physical variables strongly correlate with both global M* (Gavazzi 2009) and (Dopita 
& Ryder 1994) suggesting that both stellar mass and its surface density have significant roles in galaxy evolution. 
Strong positive correlations also exists between spatially resolved measurements of E*, and Eh2 and E24/zm- Figure 
shows such correlations (p 0.5 — 0.8) in both - Eh2 and E* - E24pm diagrams for local measurements of STIN 
galaxies. With the assumption about E* made above, we express these relations by the following power laws, 

EsFR = a-E,"\ (Al) 

EH2=b-E*", (A2) 

where (a, b) are normalization constants and (m, n) are power law indices. For each galaxy we derive these parameters 
using the OLS Bisector regression method in log space. We also obtain observed scatter (uobs) in each of these relations 
using the best-fit line. These parameters are used as inputs in the Monte Carlo simulation and are shown in Table [2j 
We use the transformation method to randomly sample stellar mass surface density from its observed distribution 
(see section 7.2 in Press et al. 1992 for a detailed account of transformation method). This method provides three 
simple steps to generate random numbers from any arbitrary distribution function: 1) construction of the normalized 
cumulative distribution function of any given variable, 2) selection of a random number from a uniform distribution 
and inspection of the cumulative function to find the match between the two, and 3) selection of the variable within 
the domain that corresponds to the match in step 2. We use this method with two additional conditions while 
sampling E* from the observed distribution. These conditions are: 1) the lower and upper limits of the samples 
surface density must be within the range of observation, and 2) the number of observed measurements and simulated 
points must be the same. With a randomly sampled E*, we generate SFR and molecular gas mass surface densities 
usin g, Y — a ■ ■ lO^(°'^obs)^ where a = {a, b}, (3 — {m, n}, and {X, Y} are the corresponding variables in equations 
|A1| and |A2| The last term introduces scatter in the relation assuming a normally distributed CTobs- Following this 
relation, we derive Y for any given X, which is then offset randomly using CTobs- We repeat this for all the points 
contributed by any given galaxy. For each galaxy we generate 5000 realizations and for each realization we compute 
stellar velocity dispersion, dynamical timescales, and sSFR. A similar Monte Carlo method has been used previously in 
two-dimensional color-magnitude fitting analysis by Dolphin et al. (2001) and in constraining the Egpp - Eh2 relation 
by Blanc et al. (2009). 

NGC 3949 has a very small dynamic range (^-^0.3 dex) in both Eh2 and E24 ^m yielding a negative correlation in 
the Eh2 - E* plane and a correlation coefficient close to zero in the Eh2 - E24pm plane. We exclude this galaxy from 
our simulation. It contributes only 2% to the ensemble of points and has virtually no impact on our results. We also 
exclude the seven outhers with E*' > 2000 Mq pc"^ from NGC 5371 data. 

The results of our simulations are presented in Figs. [7] and |8] Figure [7] highlights various observed relations 
and their reproductions from a single realization. While panels 1 and 2 demonstrates relationships between dynamical 
timescales and Eh2 and E24pm for the observed data, panels 3 and 4 show these timescales for simulated measurements. 
Simple correlation analysis strongly suggests that the simulation reproduce the observation reasonably well. Simulated 
measurements, however, show slightly larger scatter (~0.15 dex in panels lb and 3b as compared to ~^0.25 dex in panels 
2b and 4b), which is mostly due to CTobs- Since Uobs is a quadratic sum of intrinsic scatter ((Tint) and measurement error 
(cmes), in absence of (Tmes the observed scatter would always overestimate the true scatter present in any relation. 

Figure [8] shows all three major correlations for both observed and simulated measurements. For the purpose of 
demonstration the simulated data are drawn from the same realization as in Fig. [T] A comparison of the corre- 
lation coefficients of individual galaxies in Tables [T] and [2] shows that the simulation approximately reproduces the 
observed r^"' - dynamical timescale correlations. The simulation also produces r^°' - sSFR correlation comparable 
(—0.71 < p < —0.62) to the observed one {p ~ —0.64). The right panels show the normalized distribution functions of 
the correlation coefficient derived from 5000 realizations. The mean (p) and the standard deviation (cTp) of the corre- 
sponding distribution functions are (-0.08±0.02), (-0.22±0.02), and (-0.68±0.01). The mean values of the respective 
distribution functions reflect the intrinsic strengths of these correlations. For example, a small correlation coefficient 
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Figure 6. The observed S24/im - S» and Sh2 - correlations for STING galaxies. The measurements of each galaxy have different 
color code. The measurements of NGC 772 and NGC 3147 are highlighted using same color codes and legends as in Fig. [T] 



Table 2 

Input and Output Parameters of Monte Carlo Simulations 



Object 


log a 


m 


Cobs 


logb 


n 


Cobs 


N 


Corr. 


Corr. 








dex 






dex 




P 


P 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC 628 


-0.50 


0.62 


0.11 


-0.81 


0.92 


0.11 


131 


-0.23 


-0.45 


NGC 772 


-1.12 


0.85 


0.19 


0.16 


0.63 


0.19 


217 


-0.40 


-0.54 


NGC 1637 


-0.80 


0.82 


0.09 


-1.45 


1.27 


0.09 


47 


-0.72 


-0.76 


NGC 3147 


-0.20 


0.49 


0.06 


-0.65 


0.90 


0.06 


298 


-0.79 


-0.82 


NGC 3198 


-3.65 


2.13 


0.08 


-2.99 


1.89 


0.08 


18 


0.41 


-1-0.24 


NGC 3593 


-1.53 


1.08 


0.17 


-0.74 


0.89 


0.17 


141 


0.19 


-0.04 


NGC 3949 














27 






NGC 4254 


-0.14 


0.64 


0.10 


-0.58 


0.99 


0.10 


308 


-0.64 


-0.67 


NGC 4273 


-1.10 


1.12 


0.09 


-1.19 


1.22 


0.09 


103 


-0.10 


-0.16 


NGC 4536 


-0.44 


0.89 


0.15 


-0.52 


0.94 


0.15 


67 


-0.16 


-0.28 


NGC 4654 


-1.05 


0.96 


0.10 


-0.93 


1.03 


0.10 


168 


-0.20 


-0.32 


NGC 5371 


-1.41 


0.94 


0.11 


-1.50 


1.22 


0.11 


58 


-0.57 


-0.61 


NGC 5713 


-0.94 


1.07 


0.21 


-0.59 


0.99 


0.21 


220 


-0.12 


-0.20 


NGC 6951 


-1.98 


1.26 


0.18 


-1.54 


1.26 


0.18 


135 


-0.10 


-0.17 


Note. — 


Column (1): Galaxy 


name. 


Column (2), (3) I 


^ (4): Normalization constant, power- law index, 


and 


bserved scatter (in 


dex) derived 



from S24^m - relation. Column (5), (6) & (7): Normalization constant, power-law index, and observed scatter (in dex) derived from Sh2 - 
relation. Column (8): Number of independent observed measurements. Column (9) & (10): Spearman rank correlation coefficients (p) from 



and 



relations, respectively, for one realization of Monte Carlo simulation. 



(p < 0.2) derived from a large set of measurements indicates no correlation between t^™' with either Jeans timescale 
or free-fall timescale. Likewise, one can find a strong anti-correlation between T^p'and sSFR. 

The properties of the distribution functions shown in Fig. |8] depends on how scatter is introduced in Eqs. |A1| and 
|A2| While simulation of each individual galaxy takes normally distributed Cobs as an input, the observed scatter, 
however, is not symmetric and varies with surface density (see Fig. pi). An outcome of this simple choice of dobs is 
the narrow distribution functions of p and the offsets between p and observation. However, even with this sim plistic 
approach we approximately reproduce the observed strengths of various correlations using Eqs. XT and X2 which 
suggests that these relations are mostly determined by the interrelations among S24/im, Sh2, and and c*. 
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Figure 7. Observations vs. simulations for STING galaxies. The dynamical timescales as a function of Sh2i and inverse S24/xm are 
shown in panels 2 and 3. These timescales derived from a single realization of the Monte Carlo simulation arc shown in panels 4 and 5. Each 
panel shows the Spearman rank correlation coefficient (p) for the corresponding relationship. The observed and simulated measurements of 
each galaxy have the same color code as in Fig. [6] The measurements of NGC 772 and NGC 3147 are highlighted using same color codes 
and legends as in Fig. [l] 
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Figure 8. Observations vs. simulations of STING galaxies. The 
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sSFR relations for observed and 
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simulated measurements are shown, respectively, in panels (la), {2a), and (3a) on the left. The open contours (inwards) enclose 99%, 
75%, 50% and 25% of the observed distribution. The simulated measurements comes from the same realization and have the same color 
code as in Fig. [7] The right panels show the normalized distribution functions of the Spearman correlation coefficient obtained from 5000 
random realizations of the data. The mean (p) and the standard deviation {(Jp) of the corresponding distribution functions are -0.08it0.02, 
-0.22it0.02, and -0.68±0.01. The correlation coefficient for the corresponding observed relations are -0.01, -0.16, and -0.64. 



